Spin-orbital physics in the optical conductivity of quarter-filled manganites 
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Using finite-temperature diagonalization we investigate the optical conductivity o~(u>) and spin- 
orbital dynamics in the CE phase of half-doped manganites. We find o~(u) characterized by a broad 
spectrum with pronounced optical gap due to charge ordering induced by Coulomb and further 
neighbor Jahn- Teller interactions. With increasing temperature the conductivity shows a significant 
change over a wide energy range with a characteristic shift towards lower frequencies. In the low 
temperature CE phase we observe in-gap absorption due to combined orbiton-spin excitations. 

PACS numbers: 75.10.-b, 71.70.Ej, 75.30.Et, 75.40.Mg 



The large variation of physical properties of mangan- 
ites as function of doping or temperature originates from 
a complex interplay of spin, orbital and charge degrees of 
freedom, as well as the interaction with the lattice 0, • 
An important example for the control of magnetic order 
due to the orbital degree of freedom is the CE phase at 
quarter filling. Here ferromagnetic (FM) zig-zag chains 
are staggered antiferromagnetically and the occupied e g 
orbitals are oriented along FM bonds in the (a, b) planes 
[Fig-E&;b)]. Although the structure reflects the coop- 
erative action of antifcrromagnetic (AF) superexchangc 
and FM double-exchange (DE), these interactions alone 
do not guarantee its stability 0, |fj . The AF-CE_phase 
was observed in cubic (Nd,Prb /o(Sr,Ca) 1 / 2 Mn02_p] and 
in layered La 1/2 Sr 3 / 2 Mn0 4 :6|Q, LaSr 2 Mn 2 7 Qman- 
ganites. Yet CE-type orbital correlations have also been 
observed in the absence of antiferromagnetism, namely in 
the FM metallic phase of Ndi/ 2 Sr!/ 2 Mn03|2(, and also 
at smaller hole doping and are believed to be connected 
with the mechanism of colossal magnetoresistance (lfij . 
An important feature of the CE-phase is a large optical 
gap and a broad absorption maximum at 1-2 eV shifting 
towards lower frequencies with increasing temperature 
[ill H2I Il3| . Remarkably a similar spectral shift is ob- 
served in photo-excitation experiments, which reveal ul- 
trafast response times of such Mott insulating structures 

urn. 

Our aim is to investigate the electron dynamics (opti- 
cal conductivity, elementary excitations) in the CE spin- 
orbital state. It is evident that a real understanding of 
this phase is only achieved if also the excitation spec- 
tra are consistent with experiment. The strong variation 
of o~(lo) with temperature T was not addressed in previ- 
ous theoretical work [lfil Hf| . Our study employing finite 
temperature diagonalization is designed to explore also 
quantities, such as the orbital excitation spectra, that so 
far cannot be measured by experiment and investigate 
their effect on observables like optical conductivity and 
spin excitations. Our results for the optical conductivity 
reproduce for the first time the experimental trends on 
a semi-quantitative level, and predict orbital excitations 



in the charge gap of the AF-CE phase. 

Orbital degeneracy provides the option for the system 
to lower its dimensionality, which is essential for the ap- 
pearance of the CE structure 0. Our study explores 
a novel aspect of orbital physics in manganites, namely 
the role of further neighbor Jahn- Teller (JT) interactions, 
and the consequences on the dynamics. We find that or- 
bital polarization induced by this interaction is crucial for 
the appearance of CE-antiferromagnetism at low T, and 
leads to a spin response consistent with AF-CE order. 

We consider an extended DE model with charge or- 
der (CO) triggered by Coulomb repulsion and further- 
neighbor JT-interaction ||. AF order along the c axis 
and e g orbital elongation in the a/b direction confines 
the electron dynamics to (a, b) planes and thus we con- 
sider a two-dimensional (2D) model. Assuming infinite 
on-site Coulomb repulsion between e g electrons on the 
same site (U — + 00) the model has the form ||: 

<ij> CO <i> 

+ J A pJ2 S i S ! + Vj2nm i + Hoo- (1) 

(«> <y> 

The first term describes the motion of e g electrons un- 
der the constraint that each site can be occupied at most 
by one e g electron, = d\ ia {l - n-^ s ) Yl a ,{l - n^,). 
Here the index £ (a) denotes the e g orbital (spin) orthog- 
onal to £ (cr), respectively. The hopping matrix elements 
ty in the orbital basis {\x), \z)} with \x) ~ x 2 — y 2 and 
\z) ~ 3z 2 — r 2 have the standard form [l8j. The sec- 
ond term in Eq.Q stands for Hund's coupling between 
electron spin Si = X^o-a-' dl^daa'di^i and the i 2g core 
spin Sf (S c = 3/2). The third term represents the AF 
superexchange (Jaf > 0) between nearest- neighbor i 2g 
spins, while the fourth term is the inter-site Coulomb re- 
pulsion which favors checkerboard CO (V > 0). Finally 
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accounts for the cooperative JT-interaction and de- 
scribes the orbital-orbital (00) interaction between fur- 
ther Mn 3+ neighbors mediated by the Mn 4+ ion in- 
between [j| [Fig. ^c)]. Hence, (jij') denotes three neigh- 
boring Mn 3+ -Mn 4+ -Mn 3+ sites along a or b direction 
while the two-site orbital operator, Tjj/, in the {\x), \z}} 
basis has the form T ]v = T. Z TZ + 3T?T$ T VZ(T?T?, + 
T?Ty) expressed in terms of pseudospin operators: T? = 

\ E CT ("i^-n izCT ) and If = \ E CT (d\ x J- 1Z!y + d\ z Ji x ^ ■ 
This term not only leads to the CE-type orbital order but 
also induces charge alternation on neighboring Mn sites. 
The interaction Hqq originates from a JT-driven dis- 
placement of the Mn 4+ octahedra as observed in the CE- 
phase of La^Ca^MnOslljj. The nearest-neighbor 
00 interaction 2n E Ty [20( is neglected here for sim- 
plicity as it does not contribute in the case of strong CO. 

First we focus on the optical conductivity which is de- 
termined by the current-current correlation function 



cr(w) 



Nuj 
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where the cc-component of the current operator is j x = 
i e E(ij) E ?C:<T t^ S x (dl^dxv - H.c.) and 5» denotes the 
^-component of the vector connecting sites i and j. We 
evaluate this Kubo formula using the exact diagonaliza- 
tion technique for finite temperatures developed by Jaklic 
and Prelovsek 1211 for cuprates and later used also for 
manganites 0, Il8| . In order to make the problem numer- 
ically tractable, we simulated the core spins by S c = 1/2 
spins and scale Jh to preserve the splitting between high 
and low spin states. The data reported here was calcu- 
lated for Jh = 15i. Although we exploited the transla- 
tional symmetry to save storage and computation time, 
calculations were restricted to a small \/8 x \/8 cluster 
as only Sf ot = Ei subspaces (Si = Sj c + s^) can be 
treated separately, while a similar symmetry in the or- 
bital sector is absent. 

The subtle role of -ffoo is to enhance CE type or- 
bital polarization thereby supporting the AF decoupling 
of FM-chains at low T (Fig. [TJi), while for small k' DE 
dominates in all directions leading to a FM charge and 
orbital ordered state [f| . In Fig. we show the fre- 
quency and temperature dependence of emerging 
from ground states with AF-CE [Fig. I^a)] and FM-CE 
[Fig-Htb)] spin-orbital order, respectively. A characteris- 
tic feature of cr(w) is a broad absorption maximum found 
for 2t < ui < ht which shifts towards lower frequencies 
by Aui ~ 0.5i as temperature increases to T ~ 0.2i. 
Further increase of temperature leads to the melting 
of the orbital order and the closing of the charge gap 
(lu < 2t) with further shift of the high frequency broad 
structure. The shape of o-(w) and its temperature evo- 
lution is similar to the experimental data, e.g., for the 
cubic Lao.sCao.sMnOa EH and layered La-i /?Sr3/?Mn04 
compound. Assuming t = 0.4 eV ^[ one finds 
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FIG. 1: (Color online) Sketch of the AF-CE spin, orbital, 
and lattice order. In (a) large (small) spins indicate sites 
with (without) an e 9 electron, respectively. Shading indicates 
direction of DE carrier propagation. In (b) full symbols rep- 
resent occupied e g orbitals and in (c) full circles indicate the 
JT-shifted O ions, and the arrow the displacement of Mn 4+ . 




FIG. 2: (Color online) Temperature dependence of the optical 
conductivity o(uj) obtained in the phase with: (a) AF-CE 
order in the ground state, and (b) a FM-CE ground state for 
Jaf = 0.04i. Other parameters: (a) V = 0.2t, k = 0.15i; 
(b) V = t, k' = 0.05t. The spectra are generated using a 
Lorentzian broadening T — 0.2t. Inset: Mid-gap absorption 
broadened with V = 0.05t. 



the position of a broad peak structure in Fig. Reentered 
close to the experimental value 1.2 eV measured in single- 
layer Lai/2Sr 3 / 2 Mn04 A subtle test is the absolute 
scale of a(u>) which is controlled by the sum rule that 
links cr(o>) and the kinetic energy of the model 
For the maximum of the absorption (a max ~ O.I/Oq -1 ) 
we find with po = ha/e 2 E3 and the lattice constant 
a ~ 5.5A, a max pa 500 (Q.cm)^ 1 . Observed values for 
o-max are ~ 1000 (ftcm)" 1 in La 1/2 Sr3 /2 Mn04_El 
and ~ 700 (ficrra) -1 in cubic Lao.sCao.sMnOa jl3j. 

The main absorption in Fig. [5] describes excitations 
across the charge gap (A c ~ 3(V + 2k') in case of per- 
fect CO) which in the orbitally ordered state involve e g 
electron transfer mainly between bridge (ground state) 
and corner (excited state) sites within the same zig-zag 
chain. To understand the role of hopping t it is useful to 
consider an embedded three-site model with an e g elec- 
tron being either at the bridge (\B)) or at one of the 
corner (|Cl), \Cr)) sites. One finds the optical transi- 
tion at uj c ~ A c + 2t 2 /A c , between bonding ~ \B) and 
non-bonding (\Cl) — \Cr))/V2 states which have differ- 
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ent parities. With increasing T magnetic excitations are 
included in the thermodynamic average with reduced ef- 
fective hopping, whose value is dictated by the DE mech- 
anism for correlated e g electrons. Thus the reduction of 
the second order term 2t 2 /A c for these states results in a 
spectral shift towards lower u> c . Further reduction of u c 
is caused by melting of orbital and CO leading to gapless 
orbital excitations and a reduction of A c . 

In the FM spin state ct(uj) is characterized by two broad 
peaks and a well-developed charge gap [Fig. 0(b)]. Here, 
spin fluctuations are small and cr(u>) is dominated by 
charge-orbiton excitations, i.e., the pure charge excita- 
tion is well separated from two further peaks correspond- 
ing to the additional excitation of one and two orbitons. 
On the other hand, in the AF-CE spin ordered state spin 
fluctuations are substantial and inter-site charge trans- 
fer is accompanied by magnon and orbiton excitations 
leading to a broadened spectrum above the charge gap 
as well as to substantial in-gap absorption. The latter is 
completely absent in the FM state at low temperature. 

To gain a deeper insight into the structure of <t{lo) in 
the CE spin-orbital phase, we have calculated [using a 
similar finite-temperature procedure as for <r(u>)] the dy- 
namical spin, orbital, and spin-orbital response functions: 
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where A q = 5 q , T q , or C q stand for the spin, orbital, 
or combined spin-orbital operator C q = ^ k >Sk2k+ q , re- 
spectively. These response functions describe the sim- 
plest elementary excitations realized in the spin-orbital 
model. Strongly broadened collective spin excitations are 
present at low temperatures, but already at T ~ O.lt 
quasi-elastic scattering (w ~ 0) is found at all mo- 
menta (Fig.[3J). The spin spectrum at low temperatures 
is strongly anisotropic, with pronounced precursors of 
Bragg peaks present at momenta (tt, 0) and (7r/2,7r/2) 
at uj/t ~ 0.05 reflecting the spin pattern of the AF- 
CE phase [see Fig. ma)]. Here, spins S = 1 form two 
interpenetrating AF superlattices (with lattice constant 
a = 2) which leads to the main maximum in Sq(u>) at 
q = (7r/2,7r/2), while spins S = 1/2 form one AF lattice 
in the frame rotated by 45 degrees leading to the smaller 
maximum in Sq(uj) at q = (-7T, 0). Spin waves develop at 
higher energies lo ~ 0.2t due to the competition between 
AF t2 g -t2g interactions and DE mechanism. An accurate 
assessment of the spin- wave dispersion would, however, 
require the study of larger clusters. 

Orbital excitations are much more robust than the 
magnons and only weakly change with increasing tem- 
perature from O.Olt to O.lt. In the T q (o;) spectrum [see 
Fig. EJa)-(b)] one finds a pronounced peak only at mo- 
mentum q = (7r/2,7r/2) with u> ~ 1.3t at low temper- 
atures which broadens with increasing temperature and 
becomes incoherent only for T > t. Less coherent multi- 
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FIG. 3: (Color online) Dynamic spin structure factor Sq(u>) 
in the AF-CE (T = O.Olt) and paramagnetic (T = O.lt) phase, 
respectively. Parameters: V = 0.2, k! = 0.15t, Jaf = 0.04t, 
Jh = 15t. Spectra are broadened with T = 0.02t. 
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FIG. 4: (Color online) Dynamic structure factor in the {|x}± 
\z}} basis for orbital excitations Tq(o>) (a) T = O.Olt, (b) 
T = O.lt, and (c) combined spin-orbital excitations C q (u) at 
T = O.Olt. Other parameters as in Fig.0]and F = O.lt. 



peak structure is found at other momenta. The fact that 
the orbiton spectrum shows a pronounced peak only at 
(tt/2, 7t/2) momentum is a consequence of the pseudospin 
structure of the CE orbital phase [see Fig. ^b)] which 
have the same pattern as the S = 1 spin structure [see 
Fig. IHa)]. In the {|a:) ± \z)} rotated orbital basis the 
OO term [see Eq. ©] has Ising-like form, Ty ~ 3TfT7, 
and in the case of perfect Mn 3+ -Mn 4+ CO the energy of 
orbital excitation is estimated ~ 12k'. Thus, for 

k' = 0.15t one expects u> or b — l-8t while in our numer- 
ical results [see Fig. 0I a )] * ne orbiton energy is reduced 
to 1.3t due to imperfect CO. Comparing the spectrum 
of orbital [Fig. IHa)] an d spin-orbital [Fig. EI C )] excita- 
tions one finds that collective modes are formed as well. 
At energy lj ~ 0.7t and different momenta one finds a 
peak well below the orbital mode (w ~ 1.3£) which 
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FIG. 5: A FM segment showing the dominant terms in the 
AF-CE ground state with CO (a, (3 <C 1) in the limit Jh — > oo 
(left). When applying the current operator j x (right), only the 
first term would contribute in the FM case. The second term 
~ (3 in J X |GS) (describing an orbital excitation) is finite due 
to quantum fluctuations in an AF environment. Here t and 
t' (t > t') stand for two effective hoppings between corner 
and bridge sites, while full (empty) symbols stand for the 
sites occupied (unoccupied) by an e 9 electron with total spin 
(arrows) S = 1 {S = 1/2), and |0) = \S = 1, S z = 0). 



stands for a bound state |23|, |2J| while at oj ~ 1.5i a 
strong maximum is found having a character of an anti- 
bound state. This indicates that the orbital excitations, 
although much more robust than the spin ones, play an 
active role also at a lower-energy scale. At higher energies 
u> > 1.5t, however, Cq(w) and T q (u;) structure factors 
strongly differ. 

A salient feature of cr(u>) is the in-gap absorption 
(oj < A c ) in the AF-CE phase, Fig. 2(a), and the absence 
of such absorption for the FM-CE case in Fig. 2(b) at 
low temperature (T < Q.lt). We argue that the structure 
at u) ~ t is an orbiton-like collective excitation which be- 
comes infrared-active when parity is broken by AF corre- 
lations. This is illustrated schematically in Fig. O In the 
CO phase the ground state (GS) is dominated by config- 
urations with occupied bridge sites but hybridized with 
states reached by inter-site hopping (~ t) or changed by 
core spin superexchange (~ Jaf)- As we deal with d 
orbitals the hopping matrix has even parity leading to 
even GS as far as charge distribution is concerned. Con- 
sequently, the GS changed by the current operator is odd 
with respect to charge and (neglecting core spins) only 
processes across the optical gap are allowed (first line of 
ix|GS) in Fig-EJ- In the presence of AF fluctuations these 
selection rules no longer strictly hold, as spin-exchange 
leads to an admixture of spin and spin-orbital excitations 
(second and third line of j x |GS) in Fig. O which do not 
cancel by symmetry. 

This orbiton-magnon mechanism for the in-gap ab- 
sorption is supported by the fact that the intensity of 
the structure in <t(uj) at w ~ t [see inset in Fig. Ufa)] 
increases (above the background of other excitations) 
with decreasing temperature in the spin-ordered case 
(T < O.li). The momentum conservation present in op- 
tical processes can be fulfilled combining and 5* ex- 
citations with q = (±7r/2, ±7r/2) which are dominant in 
the respective dynamical structure factors [see Fig. Ofa) 
and Fig. Ufa)]. The decrease of the static structure factor 



Sq = / dujSq(ui) at q = (±7r/2, ±7r/2) with tempera- 
ture is apparently correlated with the suppression of the 
in-gap absorption in o~{u>). 

In summary, optical conductivity and spin-orbital dy- 
namics have been investigated in the CE phase of charge- 
ordered manganites at quarter filling. Two contributions 
in our theory determine the charge gap: (i) elastic forces 
and (ii) Coulomb interaction. This is important to ex- 
plain the recently discovered pico-second switching of op- 
tical spectra, which can be understood by the instanta- 
neous reduction of the Coulomb contribution due to the 
suppression of CO by the ultra-short pulse, as well as 
the fast recovery after the pulse. We have shown that 
the cooperative JT coupling between further Mn neigh- 
bors, stabilizing the CE-type correlations in our model, 
can explain semi-quantitatively the optical data and its 
temperature dependence in half-doped manganites. Spin 
and orbital excitations lead to collective modes active in 
optical experiments, and are responsible for characteris- 
tic differences between AF- and FM CE phases. 
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